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ABSTRACT 

We show that different stellar-mass estimation methods yield overall mass scales that disagree by 
factors up to ~2 for the z = galaxy population, and more importantly, relative mass scales that 
sometimes disagree by factors >3 between distinct classes of galaxies (spiral/irregular types, classical 
E/SOs, and E/SOs whose colors reflect recent star formation). This comparison considers stellar mass 
estimates based on (a) two different calibrations of the correlation between if -band mass-to-light 
ratio and B — R color (Bell et al., Portinari et al.) and (b) detailed fitting of UBRJHK photometry 
and optical spectrophotometry using two different population synthesis models (Bruzual-Charlot, 
Maraston), with the same initial mass function in all cases. We also compare stellar +gas masses 
with dynamical masses. This analysis offers only weak arguments for preferring a particular stellar- 
mass estimation method, given the plausibility of real variations in dynamical properties and dark 
matter content. These results help to calibrate the systematic uncertainties inherent in mass-based 
evolutionary studies of galaxies, including comparisons of low and high redshift galaxies. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 

The modern trend toward studying galaxy properties 
as a function of mass rather than luminosity has led 
to remarkable advances in our understanding of galaxy 
evolution, making the calibration of mass estimation 
techniques a high priority. Based on comparisons of 
dynamical and ste llar- mass (Md vn and M» ) estimates, 
iDrorv et all (|2004f) and lRettura et all (|2006f) argue that 
multi-band photometry alone can provide accurate M, 
estimates that correlate well with M& yn . Even better, 
modeling of the correlation between optical colors and 
stellar mass-to- light ratios M*/L suggests that factor-of- 
two accuracy in M* may be achievable with just three 
filters, especially w hen using optical colors to infer an I- 



' or if -band M*/L dBell fc de Jondl200l 



2004). However, recent work (Maraston 



Portinari et al.l 
2001 " cautions 



that M* estimation may be more complicated than pre- 
viously assumed, as young stellar populations may con- 
tribute substantially to not only optical but also near- 
infrared light, via thermally pulsing asymptotic giant 
br anch (TP-AGB) sta r s. Us in g these models and those 
of iBruzual fc Chariot! IpOOlh . Ivan der WeleTall (|2006h 
find substantial inconsistencies between Md yn and -M* 
estimates at low and high z that appear only when near- 
IR photometry is used. 

To date, empirical examinations of these issues have 
relied on mixed data sets, making it hard to isolate sys- 
tematics in M* estimation from evolution between low 
and high redshift galaxies and/or effects of inhomoge- 
neous data. Here, we take advantage of the high-quality, 
uniform dat a available for the N earby Field Galaxy Sur- 
vey fNFGS. Dansen et al.ll200 0b), including photometry, 
spectrophotometry, and gas and stellar kinematics, to 
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evaluate M* estimation techniques. Our sample allows us 
to explore effects of stellar population age on M* estima- 
tion at a single redshift, as it includes late-type galaxies, 
classical E/SO galaxies that fall on the red color- M* se- 
quence, and galaxies with E/SO morphologies that fall on 
the blue color- M* seque nce due to recent star for mation 
( "blue-sequence E/SOs" iKannappan et alll2006al here- 
after KGB). 

2. METHODS 

The NFGS provides a broadly representative galaxy 
sample spanning a wide range of luminosities and mor- 
phologies. For M* estimation, we analyze 141 NFGS 
galaxies with UBRJHK photometry and in tegrated (slit- 
scanned ) optical spectrophotometry from iJansen et al.l 
(|2000bllal ) and the Two Micron All-Sky Survey Ex- 
tended Source Catalog (2MASS XSC, I Jarrett et ail 
2000); see KGB for sample selection details. We 
also check our results using ugr photometry for 92 of 
these galax ies, taken from the Sloan Digit al Sky Sur- 
vey (SDSS, |AdeIminlMcCirt^irin [2006|). Photome- 
try a nd spectra are correc ted for foreground extinction 
using Schlegel et al.l (Il998f) and the Galactic extinction 
curve of lO'Donnell (|1994D . We do not apply internal 
extinction corrections to the data used for mass deter- 
mination (though corrections based on the method of 
iTullv et al. 1998 are used incidentally for defining the 
red and blue sequences, with no effect on our mass er- 
ror budget; see KGB). However, dust is either included 
in our modeling or, in the case of color- relations, 
neglected following standard practice. We add 0.1 mag 
in quadrature to the catalogued internal magnitude un- 
certainties for all passbands to account for systematic 
uncertainties in foreground exti nction correction s, auto- 
mated 2MASS photometry fsee iBell et "all 120031 . here- 
after B03), and relative photometric zero points. For 
the spectra, we add relative flux calibration uncertain- 
ties (typically 6%, but up to 9% outside 4000-6800 A; 
R. Jansen, priv. comm.) in quadrature to the formal 
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uncertainties. 

Our "reference" M* values are computed by fitting the 
photometry and spectra to a discrete grid of stellar popu - 
lation synthesis models from lBruzual fe Chari ot ( 2003), 
scaled to a "diet Salpeter" IMF as used by B03. We 
combine two simple stellar populations (SSPs) in vary- 
ing mass fractions (100:0%, 90:10%, 80:20%, etc.), where 
we have normalized each SSP to M* = 1 Mq. Each 
individual SSP has one of eight ages (0.025, 0.1, 0.29, 
0.64, 1, 2.5, 5, 11 Gyr) and three metallicities (0.4Z Q , 
Z Q , 2.5Z q ), and each combination of SSPs has one of 
eleven dust optical depths (Tv, ga s = 0, 0.12, 0.24,..., 1.2). 
We derive model photometry by convolving NFGS (stan- 
dard Johnson-Cousins), SDSS, and 2MASS filter profiles 
with Bruzual -Charlot model spectra and adding attenu- 
ation using a iCalzettJ ([2001) law. The code scales each 
model to the observed photometry in L© units, yield- 
ing the estimated M* for that model, then computes 
likelihoods oc e _x »»=™«/ 2 for the entire grid of models. 
In a first pass, we fit only the photometry, redshifting 
the models to match the individual galaxy spectroscopic 
rcdshifts and comparing with redshift-zero models to de- 
termine fc-corrections. In a second pass, the likelihood- 
weighted average fc-corrections are applied to the input 
photometry, and the code fits both the photometry and 
the de-redshifted spectra to a fixed set of models in the 
rest frame. We mask emission lines, limit the spectral 
range to 3800-7000 A, and convolve the model spectra 
to the 6A resolution of the NFGS spectra. As the spec- 
tra lack absolute flux calibration, their scale factors are 
allowed to vary freely. The likelihood of each model is 
the product of the likelihoods inferred from the photom- 
etry and the spectra, so the y 2 terms sum in the expo- 
nent and can be weighted to set the relative influence 
of the spectra and photometry. We adopt Xoveraii = 

max (X 2 pho V n dof) + Xlpec-raw/ 1000 + Xspec-norm/ 1000 > 

where the latter two terms are contributions from fits to 
the raw and continuum-normalized spectra and ridof is 
the number of degrees of freedom in Xphov nOTma lly nve 
when fitting six filters (losing one to the scale-factor de- 
termination) . Once the photometric data are reasonably 
well fit (Xphot — n dof), the likelihoods are affected only 
by the spectra. Otherwise, the likelihoods are equally in- 
fluenced by the reduced-^ 2 values of the spectra and pho- 
tometry, because the ratio of the number of data points 
is ^500 (where the strong covariance between X% ec -raw 
and xipec-norm justifies treating them as a joint x1 pe c 
term). Following iBundv et alj (|2005i ), we adopt the me- 
dian of the likelihood distribution binned over logM* 
rather than the best fit to determine the final M*, and 
wc estimate uncertainties from the 68% confidence inter- 
val in logM* (binning in 0.02 dex intervals). 

M* estimates based on the B — R vs. M„/Lk re- 
lation are derived fr om the calibrations of B03 and 
iPortinari et~aT1 (|2004 hereafter P04). The B03 cali- 
bration is based on a global linear fit to M*/Lk vs. syn- 
thetic B — R for a large sample of galaxies with ugrizK 
data, where each galax y is fitted with PEGASE popula- 
tion synthesis models ([Fioc fc Rocca-Volmerangg|1997l) 
to find the best-fit metallicity and exponential star for- 
mation history (SFH), which may be decaying, constant, 
or rising. The P04 calibration is predicted from chemo- 



photometric models of galactic disks, which include TP- 
AGB stars. We use P04's Salpeter IMF calibration, mul- 
tiplying the resulting masses by 0.7 to match the diet 
Salpeter IMF scale of B03. Technically, P04 limit their 
calibration to B — R = 0.95-1.45, and our use of the 
relation sometimes extends outside this range. Factor 
of two uncertainties are predicted for color-based mass 
estimation, primarily due to variations in SFH but also 
due to the neglect of internal reddening and extinction 
(expected t o vary mainly along t he B — R vs. M^/Lk 
relation, see lBell fc de Jongll2001l) . Because color- M^/L 
relations do not provide a self-consistent way to compute 
fc-corrections, we fc-correct the input magnitudes using 
our standard method described above. 

When comparing and Md yn , we apply more restric- 
tive sample selection criteria. After rejecting galaxies 
flagge d as morphologically peculiar by iKannappan et afl 
(2002, hereafter KFF), we define two subsamples for 
which mass estimates should be robust: (1) 38 spi- 
ral/irregular galaxies with both HI data from the H yper- 
Leda homogenized Hi catalog ([Paturel et al.|[2003l ) and 
optical emission-line rotation curves passing the quality 
criteria of KFF, with the latter also having asymmetry 
<10% and extent >l-3x the B-band half-light radius rf 
(jKannappan fc Bartonll2004D : and (2) 26 E/S0 galaxies 
with optical (Mg-triplet regio n) stellar velocity disper- 
sions i n the NFGS data b ase (jKannappan fc Fabricantl 
1200 It IKannappan et al.l l2006bl ). initially measured 
within r^/4 using th e Four ier-space fitting code of 
Ivan der Marel fc Franxl (|1993l ) and rescaled to the R- 
band half-light radius rf- using eqn. 1 of lCappellari et all 
(2006). We require the rescaled dispersion to satisfy 
<r r R > 120 km s _1 to ensure negligible rotation correc- 
tions (e.g., Fig. 8 of KGB). 

For late types, gas masses are computed from the HI 
line flux with a helium-mass correction factor of 1.4 and a 
type- and mass-dependen t molecular-gas co rrection fac- 
tor of 1.06-1.4 (based on lCasoli et al.lll998l ). We adopt 
an uncertainty of 50% in the total gas mass. For the 
E/SOs, which are all massive given our cut in o t r , we 
assume gas masses are negligible (the largest measured 
gas-to-stellar mass ratio in this subsample is < 15%). 

Dynamical masses are computed for the E/S0 sub- 
sample using Mdyn = 5r^a^ R) as recommended by 
iCappellari et al.1 ((2006). We add 5% in quadrature 
to the uncertainties in for profile extrapolation er- 
rors, and 5% and 7%, respectively, to the uncertain- 
ties in a r n for template mismatch and aperture rescal- 
ing errors. For the spiral/irregular subsample we use 
Mdvn = 2.5r^K 2 ^, assuming V rot scales like \f2a (e.g., 
iBurstein et al.lll997f ). Here V ro t is 0.5x the inclination- 
corrected Wy linewidth parameter of KFF. 

We assume H = 70 km s _1 Mpc -1 and d — cz/H . 
Note that our error bars do not include distance errors 
because distances are largely irrelevant to comparisons 
between mass estimates. 

3. RESULTS 

Various estimation methods are compared in 

Fig. 1, with the ratio between the alternate and refer- 
ence 's plotted as a function of reference M* in each 
panel. Fig. 2 compares Af„( +g ) estimates against Md yn 
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Fig. 1. — Comparison of M» estimates obtained by various methods, with symbols coded to show spiral/irregular types (crosses), 
red-sequence E/SOs (gray squares), and blue-sequence E/SOs (black dots). Each panel shows the ratios between M*'s computed by an 
alternate method, as noted, and our reference M*'s computed by fitting Bruzual-Charlot models to UBR JHif+spectrophotometry, with 
ratios plotted as a function of reference M» . 



in the same format, where the notation M*( +9 ) indicates 
a gas mass correction for late types (§ HJ. We stress 
that differences between methods do not imply that one 
or the other is correct, and differences relative to Md yn 
must be interpreted carefully, given the potential for real 
variation in dark-matter content or structural properties. 

Figs, fa-c and 2a-c test the robustness of our refer- 
ence M*'s against substitution of ugr for UBR, omis- 
sion of the spectra from the fits, and exclusion of the 
near-IR JHK data. Replacing UBR with ugr, a slight 
offset appears in Fig. fa. NFGS photometry appears 
more reliable: Fig. 2b demonstrates that M„( +s ) esti- 
mates based on SDSS photometry yield greater scatter 
relative to Md yn for late-type galaxies, and inspection of 
the fits reveals that SDSS data have a fairly high rate of 
catastrophic errors relative to 2MASS and NFGS data, 
perhaps due to systematic errors in defining galaxy aper- 
tures or profiles. Note that the apparent improvement in 
scatter for E/S0 galaxies in Fig. 2b is probably fortu- 
itous: the open squares mark galaxies from Fig. 2a that 
do not have SDSS data, showing that they tend to be the 
galaxies with the largest scatter. M*'s obtained with and 
without spectra are closely consistent (Fig. f b) and com- 
pare similarly with dynamical mass (not shown) . Taking 
advantage of this result, we have verified that our refer- 
ence M*'s are robust to using a finer resolution in mass 
ratio between the two SSPs: 100:0%, 98:2%, 96:4%, etc, 
where for computational efficiency only the photometry 
is fitted. Finally, omission of near-IR data (Fig. Ic) pro- 
duces generally consistent results, with a few outliers. 
The outliers are all late-type galaxies with fairly low 
surface brightness, whose 2MASS magnitudes may be 
underestimated. Alternatively, in such bursty systems, 
M*/L may be overestimated without the IR data to an- 
chor the fits. These systems do not show large shifts from 



Fig. 2a to 2c because their M*( + q ) is gas-dominated. 

Substituting iMarastonl |2005) models for Bruzual- 
Charlot models, stronger differences emerge (Figs, fd-f 
and 2d-f). Fig. Id shows a factor of two difference in the 
relative M* scales of high-mass red-sequence E/SOs com- 
pared to both blue-sequence E/SOs and late-type galax- 
ies. With spectra omitted (Fig. fe), this difference grows 
to a factor of three, whereas with near-IR JHK data omit- 
ted (Fig. If), it disappears entirely, though a small overall 
scale difference remains (with Maraston models yielding 
<f.3x lower M*). Excluding near-IR data, M*( +g ) es- 
timates based on Bruzual-Charlot and Maraston models 
compare nearly identically to M^ yn (Fig. 2c-d). The 
differences when near-IR data are included come almost 
entirely from the Maraston models and cause shifts rel- 
ative to Mdyn in Fig. 2e-f. However, these shifts are 
within the uncertainties and might be physical: in a 
hierarchical scenario, late-type galaxies may have both 
more dark matter and more scatter in dark matter con- 
tent than early- type galaxies. Also, blue-sequence E/SOs 
may have basic structural differences from red-sequence 
E/SOs that lead to overestimated Md yn , if their a and/or 
r values are elevated due to incomplete post-merger evo- 
lution or disk-building processes (KGB). 

Figs, f g-i and 2g examine M* estimates based on color- 
M*/L relations taken from B03 and P04. Using the B03 
relation as given, we find strong disagreement compared 
to our reference M*'s in both overall mass scale and rel- 
ative scales between different galaxy classes, by factors 
>3. However, inspection of the data from which the B03 
relation was determined (their Fig. 20) reveals that their 
linear fit is skewed by a two-component distribution, con- 
sisting of a dominant linear locus and a cloud of outliers 
with blue colors and high M*/L. These outliers may be 
analogous to the outliers we find in Fig. fc & f or may 
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Fig. 2. — Comparison of M*(_|_ g ) (gas-mass corrected for late 
types only, § [2} and Md yn for various methods of M* estimation. 
Symbols are as in Fig. 1, with the addition of an open square 
symbol in panel 6 to mark E/SOs from panel a for which SDSS ugr 
data are not available. Least-squares fits and typical error bars are 
shown for early and late types, with blue-sequence E/SOs included 
with red-sequence E/SOs. 

reflect low metallicities (B03). In any case, the overall 
trends can be harmonized if we refit B03's color- M*/L 
relation using only the primary linear locus. Fig. lh 
demonstrates good agreement with only a small over- 
all scale difference when we adopt the modified relation 
logM»/L 7< = -0.616 + 0.34(5 - R) for (B - R) > 1.2 
andlogM„/L A - = -0.808 +0.5(B-R) for (B-R) < 1.2, 



fitted by eye to their Fig. 20. These formulae include a 
factor of 1.2 to convert between the M* scales of the PE- 
GASE and Bruzual-Charlot models; this offset is noted 
by B 03 when compa r ing th eir fit to earlier predictions 
from iBell fc de Jongl (|2001l) and can be estimated from 
their Fig. 20. In addition to adopting a modified color- 
M*/L relation, Fig. lh adjusts the color-based M*'s for 
the influence of dust and starbursts. B03 estimate that 
dust and starbursts, if modeled, would add ~15% and 
~10% respectively to their M* estimates, so we boost 
M*'s for galaxies with detected Ha by 15% and M*'s 
for all galaxies except massive E/SOs (M* > 10 11 M Q ) 
by a separate 10%. With this matched-B03 calibration, 
we find good agreement between the two color-based 
M* estimation methods (Fig. lh-i), with overall scales 
~1.5-1.8x higher than our reference M*'s, possibly due 
to differences in assumed SFHs and/or photometric zero 
points. There may be a slight tendency for the P04 cal- 
ibration to give higher M* to younger galaxies, which 
reduces the offset between E/SOs and late- type galaxies 
in comparison to Md yn (Fig. 2g). This apparent improve- 
ment in the match between M*(+ g ) and Md yn should be 
taken with a grain of salt, as it reflects the loss of in- 
formation on why galaxies are blue or red (age, dust, 
metallicity) . 

In summary, our results demonstrate systematic un- 
certainties in M* estimation corresponding to factors up 
to ^2 overall and >3 between distinct galaxy classes, 
even using our modified B03 calibration and matched 
IMFs. Outliers affect B03's original color- M*/Lk cali- 
bration and also emerge when we compare results from 
stellar population modeling with and without IR data 
(note both B03 and this work rely on 2MASS). More 
generally, M* estimates are highly sensitive to IR and 
spectral information when using Maraston models, es- 
pecially for red- and blue-sequence E/SOs. An Occam's 
razor argument might justify preferring Bruzual-Charlot 
models, which yield consistent results with or without 
spectra or IR data. However, we are unable to find a 
strong physical argument for preferring a particular set 
of models based on comparisons with M ( i yn , and we cau- 
tion that agreement between M*( +ff ) and Md yn is not 
by itself proof of better M* estimation, given evidence 
for variations in dark matter content, dynamical state, 
and age/dust /metallicity. We conclude that mass-based 
evolutionary studies of galaxies should explicitly consider 
the potential effects of systematic errors in M*, partic- 
ularly when analyzing young and old galaxies together 
across galaxy classes or between low and high z. 
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